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ABSTRACT 

A Doppler image of the ZAMS late-type rapidly rotating star LO Pegasi, 
based on spectra acquired between 12 and 15 September 2003, is presented. 
The Least Square Deconvolution technique is applied to enhance the signal-to- 
noise ratio of the mean rotational broadened line profiles extracted from the 
observed spectra. In the present application, a unbroadened spectrum is used 
as a reference, instead of a simple line list, to improve the deconvolution tech- 
nique applied to extract the mean profiles. The reconstructed image is similar 
to those previously obtained from observations taken in 1993 and 1998, and 
shows that LO Peg photospheric activity is dominated by high-latitude spots 
with a non-uniform polar cap. The latter seems to be a persistent feature as it 
has been observed since 1993 with little modifications. Small spots, observed 
between ~ 10° and ~ 60° of latitude, appears to be different with respect to 
those present in the 1993 and 1998 maps. 

Key words: stars: activity - stars: atmospheres - stars: late-type - stars: 
magnetic fields - stars: spots - stars: individual: LO Pegasi. 
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Instituto de Astrofisica de Canarias. 
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1 INTRODUCTION 



Late-type rapidly rotating stars show phenomena characteristic of solar-like magnetic ac- 
tivity, such as photospheric spots, chromospheric line emissions, coronal X-ray and radio 
emissions, as well as flaring activity. Their study has become important to improve our 
understanding of the effects of magnetic fields in stellar atmospheres, to provide further 
constraints on dynamo models of magnetic field generation, and to understand magnetic 
braking of stellar rotation. 

LO Peg (BD+22°4409) is a chromospherically active star of spectral class K5V-K7V, 
detected during the ROSAT WFC E UV all-sky su r vey an d the Extreme Ultraviolet Explorer 



(EUVE) survey. The first study by 



Jeffries et al. 



(119941 ) showed a high lithium abundance 



([Li/H]=1.30±0.25) and a Ha line in emission. They concluded that the star is young (with 
an age of approximately 30 Myr) and a likely member of the Local Association. 

LO Peg is an interesting target to study stellar activity because it is among the latest 



Lister et al. 


1999: 


Barnes et al. 


2005) 



specifically addressed the imaging of its surface (i.e., 

For stars as faint as LO Peg {my = 9.2), the achievement of a sufficient signal-to- 
noise ratio for Doppler imaging represents the most serious problem because the exposure 
time must be kept within a few tens of minutes to avoid the smearing of the reconstructed 
map due to fast stellar rotation. Therefore, we apply the method of LSD (Least Squar e 



Deconvolution) to increase the signal-to-noise ratio of LO Peg spectra (IDonati et al.lll997l ). 



This method makes use of a large number (up to ~ 1500 — 2000) of photospheric lines 
contained in an echelle spectrum to extract a mean line profile with a signal-to-noise ratio 
sufficiently high to apply Doppler Imaging techniques. In such a way, it is possible to broaden 
the sample of Doppler imaging candidates by making accessible a greater volume of space 
including not only faint field objects, such as LO Peg, but also lower main-sequence stars in 
nearby open clusters. 

Our specific implementation of the LSD approach is describe d in Sect. I3.1land introduces 



some improvements with 



(119971 1 along the lines of Barnes! (120041 ) . 



respect 



to th e original procedure by 



Donati fc Collier Cameron 



Lister et al. 



(Il999l ) based on obser- 



The first Doppler images of LO Peg were obtained by 
vations taken in August 1993. They found two main regions of spot coverage: a high latitude 
spot or polar crown with a starspot concentration towards longitude 70°, and a low-latitude 
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belt centred at a latitude of 25° ± 10°, with very lit tl e spo t coverage in the mid-latitudes 



in between. Later images obtained by 



Barnes et al. 



(120051 ) . based on data of July 1998, 



showed again a polar spot with a somewhat weaker longitude non-uniformity and several 
appendages extending from the polar cap down to a latitude of ~ 15°. Remarkably, neither a 
low-latitude band of spots nor an intermediate band free of spots were f ound. The differences 
in the spot patterns found in the two investigations were explained by iBarnes et al.l (120051 ) 



as a consequence of th e higher value of the st ellar v sin % adopted in the reconstruction made 



by 



Lister et al. 



(1l999l ). As a matter of fact, 



Barnes et al 



2005 ) reco nstructed a new map 



Lister et al. 



(Il999l ) and showed that the 



with their value of v sin i using the line profiles by 
resulting 1993 images were very similar to those obtained for the 1998 season, i.e., showing 
only high-latitude spots in addition to the polar cap. This result shows once again that the 
estimate of stellar parameters is a crucial task to obtain reliable Doppler maps of stars. In 
view of their a ccuracy and to warr ant a better comparison with previous maps, we adopt the 



parameters of 



Barnes et al. 



(120051 ). Our spectra were acquired in September 2003, therefore 
they allow us to perform a relatively long-term study of the LO Peg spot pattern in combi- 
nation with those previous studies. In particular, the prese nce of high-latitude spots appear s 



to be a crucial feature to test stellar dynamo models (e.g., 



Bushby 



2003; 



Covasetal 



20051 ). 



The surface differential rotation of LO Peg was determined by IBarnes et al.l (120051 ) by 
tracing the shear motion of starspots along their sequence of observations that extended 
for seven days. They found an equatorial acceleration with an equator-pole lap time of 
181 ±35 days, not too different from the Sun that has an equator-pole lap time of 120 days. 
However, the relative amplitude of the surface differential rotation is about two orders of 
magnitude smaller than in the Sun, i.e., ~ 0.002, given the much shorter rotation period 
of LO Peg. Unfortunately, our sequence of observations spans only four nights, thus an 
independent determination of the surface shear is not possible. On the other hand, the 
limited time interval of our observations with respect to the estimated equator-pole lap time 
makes differential rotation unrelevant for our analysis. 
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Table 1. Journal of the observations of LO Peg. 



Date 


UT start 


UT end 


Exp. time (s) 


No. frames 


S/N 


12 September 2003 


22:35 


00:13 


1800 


4 


50-80 


13 September 2003 


20:55 


04:06 


1800 


13 


50-80 


14 September 2003 


20:23 


04:00 


1800 


12 


50-80 


15 September 2003 


20:02 


04:18 


1800 


15 


50-80 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

LO Peg was observed on the nights between 12 and 15 September 2003 with the high- 
resolution spectrograph SARG0 at the 3.58-m Telescopio Nazionale Galileo, located at the 
Roque de Los Muchachos Observatory. SARG uses a R4 echelle grating (31.6 groves/mm) 
in a quasi-Littrow mode and has a dioptric camera to image the cross-dispersed spectrum 
onto a mosaic of two 2048 x 4096 EEV CCDs. 

We collected 44 spectra in total, obtained with the yellow grism cross disperser (with a 
300 grooves/mm echelle grating and a blaze wavelength of 589 nm). Each spectrum consists 
of 54 extractable orders with a wavelength coverage ranging from 462 to 792 nm and has a 
resolution of 86, 000. Table[T]shows the journal of the observations for the four nights. During 
the third night, one spectra was taken with a shorter time exposure, and after reduction 
it showed a too low signal-to- noise ratio and was discarded. For each night, flat fields and 
Ar-Th calibration lamps were also observed. 



We did not observe any template star, contrary to iLister et al.l (119991 ) who observed, 
with the same telescope and spectrograph used for LO Peg, the K5 dwarf Gl 673 (having 
an effective temperature similar to LO Peg, but a significantly smaller wsini) and the Ml 
dwarf Gl 649 (with a temperature similar to that expected for the starspots on LO Peg, 
i.e., T e g ~ 3500 K). For the reconstruction of the LSD line profiles, we adopted the LSD 
profile of Gl 673 as a template extracting its spectrum from the Elodie archivcl (for details 
see Section l3~2i) . 

Each spectrum of LO Peg was obtained with an exposure time of 1800 s, giving a signal- 
to- noise ratio ranging from 50 to 80, at the cost of some phase smearing since the star 



1 For details on SARG, see: http://www.tng.iac.es/instruments/sarg/ 



2 



See the online archive at http://atlas.obs-hp.fr/elodie/ 
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rotates by ~ 18° during each exposure. This effect actually sets the surface resolution of our 
Doppler maps, so that the use of the Elodie template, that has a resolution of 42, 000, does 
not introduce any further degradation of the map resolution. 



2.2 Data Reduction 



Data reduction was performed using IRAFg (Image Reduction and Analysis Facility), sepa- 
rately for the blue and red CCDs. Rejection of hot and b ad pixels is particularly importan t 
for our analysis as they can affect the LSD mean profiles (IDonati fc Collier Cameron! 1 199 71 ). 
Through a careful inspection, they were identified and replaced by performing a bi-linear 
interpolation along lines and columns using the nearest neighbour good pixels. We also dis- 
carded the wavelength range at the red end of each spectrum (above 6800 A) , badly affected 
by telluric lines. Pixel-to-pixel gain variations were removed using flat-field lamp exposures. 
Images were then corrected for scattered light. 

The echelle orders were extracted defining apertures on the best exposed LO Peg frames. 
All the other spectra were extracted using the same frame apertures. Wavelength calibration 
was finally performed using Th-Ar lamp spectra. 

The shape of the continuum in the echelle orders was fitted by a 3 rd -order spline func- 
tion. To guide us in the selection of the continuum intervals for the normalisation, we broad- 
ened a synthetic sp ectrum with T eff = 4500 K, log g = 4.5 (cm s -2 ), and [Fe/H] = 0.5 



( iCoelho et al.ll2005t) with a rotational profile corresponding to the LO Peg v sin 2 of 65.84 
km s _1 ( iBarnes et al.l 120051 ) and iteratively compared it with the continuum-normalised 
spectrum of LO Peg until a satisfactory agreement was obtained. 



For phasing our data, we used the ephemeris of 
HJD of maxima = 2449909.8059 + 0.4236 x E. 



Robb fc Cardinal! fll995h : 



(1) 



3 DATA ANALYSIS 

3.1 Least Square Deconvolution method 

Doppler imaging of rapidly rotating stars (v sini > 60 — 70 km s _1 ) is hampered by the 
shallow depths of the line profiles due to rotational broadening. This makes it difficult to 



3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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measure the distortion induced by photospheric spots and trace their motion along the line 
profiles as the star rotates. A very high signal-to-noise ratio 500 — 700) is therefore needed 
to perform Doppler imaging on such objects and this is difficult to obtain, especially for a 
star as faint as LO Peg, because the exposure time is limited by the requirement that the 
rotation of the star during the acquisition of a spectrum be smaller than the resolution in 
longitude, in order to avoid a blurring of the resulting map. 



To solve this kind of problems, 



Donatietal 



(119971 ) and bonati k Collier Cameronl (119971 ) 



introduced the Least Squares Deconvolution method (hereinafter LSD) that allows to ex- 
tract a mean line profile with an enhanced signal-to-noise ratio from an entire spectrum by 
means of a suitable deconvolution technique. Its basic assumption is that the o bserved spec 
trum is giv en by the convolution of an unbroadened spectrum, represented in iDonati et al 
(119971 ) and IDonati fc Collier Cameronl (119971 ) by a series of Dirac delta functions, with a 
rotational broadening profile, assumed to have the same shape for all spectral lines. More- 
over, limb-darkening is assumed to be the same for all the spectral lines, i.e., independent of 
wavelength. These simplifying assumptions are not suitable for spectral lines with an equiv- 
alent width larger than, say, ~ 0.25 A, because their profiles are significantly affected by the 
chromospheric temperature increase and cannot be assumed to be similar to those of weaker 
lines. Therefore, the LSD technique can be applied only after removing spectral intervals 
containing chromospherically sensitive and strong spectral lines (see below for details). 

In the present work we used an improved version of the LSD technique, which substi- 
tutes a synthetic spectrum to a line-list in the deconvolution procedure. Using a synthetic 
spectrum instead of approximating the unbroadened spectrum as a sequence of Dirac delta 
functions, allows to take into account the finite widths of the local line profiles resulting 
from, e.g., therm al Dopp l er bro adening and microturbulence. Our approach is the same as 
that proposed by iBarnesi (12004 ) and we refer the reader to that paper for more details. 

The rotational broadening profile contains the distortions induced by the surface bright- 
ness inhomogeneities. By convolution with the unbroadened spectrum, such distortions are 
reproduced along the profile of each spectral line. In other words, the basic assumption of 
the technique is that all the spectral lines repeat the same information on the rotational 
broadened profile. Therefore, it can be extracted by a suitable deconvolution, when the un- 
broadened spectrum is known. The advantage of such an approach is that individual line 
profiles are combined to enhance the signal-to-noise ratio of the extracted rotational broad- 



Doppler imaging of LO Peg in 2003 7 



ened profile, thus allowing us to perform Doppler imaging by using the sequence of such 
profiles. 

The mathematical formulation of the modified LSD method used in this paper is as 
follows. The flux in the spectrum at the wavelength A, .F(A), normalized to the value of the 
corresponding continuum F C (X), can be written as: 



where S(X) is the normalized spectrum of a non-rotating star with the same atmospheric 
parameters, P ro t the rotational broadened profile, v r the radial velocity, ranging from — v sin % 
to v sini, with v the equatorial velocity of rotation of the star, i the inclination of its 
rotation axis, and c the speed of light. The rotational broadened profile is normalized so 
that: flvsfai Piat{vT)dv T = 1. Equation (j2J) states that the observed normalized flux J-(X) is 
obtained by integrating the local, normalized specific intensity S(X), Doppler shifted at the 
radial velocity of the considered point, over the stellar disk. Since the star is assumed to 
rotate rigidly, the surface integration can be transformed into an integration over the radial 
velocity. 

If the brightness distribution over the photosphere of the star is not homogeneous, specif- 
ically if dark spots are present, the rotational broadened profile can be written as: 



where P u is the profile of the unperturbed photosphere and P s that of the spots, and the 
weights w u and w s are functions of the rotation phase that depend on i, vsini, the limb- 
darkening coefficient, the area and co-ordinates of the spots, and the ratio of their specific 
intensity to that of the unperturbed photosphere. To simplify our analysis, we assume that 
the spot line profile is simply a scaled version of the unperturbed line profile, i.e., P s = C S P U , 
where C s is a scaling constant whose value will be adjusted to warrant the best fit to the 
LSD line profiles. 

To transform the integration appearing in Equation (j2J) into a discrete summation for 
computation purposes, it is useful to consider the values of the local spectrum 5(A) along a 
discrete set of wavelength values, uniformly spaced in radial velocity, i.e.: Si = S(Xi), where 
Aj = (l + Aj_i, and At> r is the radial velocity interval the value of which is taken to 
sample adequately the unbroadened line profiles. In the present analysis, we choose Av r = 0.5 
km s _1 , i.e., significantly higher than the resolution of the observed spectra because we use 
a unbroadened synthetic spectrum computed at a resolution of 200, 000 by means of an LTE 




(2) 




(3) 
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atmospheric model with T e g = 4 500 K, log q = 4.5 



cm s 2 ), microturbulence of 1.0 km s 1 



and abundance [Fe/H]= 0.5 (see 



Coelho et al. 



2005|). 



We recast Equation (J2J), expressing the integration with the extended trapezoidal rule: 



N 



(4) 



Fi — X! a kSi-kProtkAv T , 
k=0 

where the index k specifies the radial velocity according to: v r k = — v sin 2 + Av r x k, 
N = 2 -^ 11 , with N + 1 being the number of points along the rotational broadened profile, 
and the a^'s are the numerical coefficients appearing in the extended trapezoidal rule: 
\ for k = 0, 

1 for < k ^ -/V - 1, (5) 
N. 

The rotational broadened profile P rot k is specified by a set of N + 1 values, whereas the 
unbroadened spectrum consists of a much larger number of points M, as numbered by the 
index i. Therefore, Equation (J3J) gives a system of M linear equations in iV+1 unknowns, with 
M>iV + l, that is to be solved by the method of the least squares. Specifically, the solution 
can be comp uted by means of t 



| for k 



Press et al 



as detailed in 
can be expressed as 

O'kSi—k 



A 



ik 



ie Singular Value Decomposition method (hereinafter SVD), 
(jl992l ). The elements of the design matrix for the problem at hand 

(6) 



and the vector of the data as: 
7 F< 

h = —, 



(7) 



where <7j is the standard deviation of the observed normalized flux Ti at wavelength \. 
In addition to the least square solution of the system (BJ, the SVD method provides the 
covariance matrix among the N + 1 values specifying the rotational broadened profile P ro t-,k, 
that allows us to derive the errors on the mean profile and its effective resolution. This is 
particularly important in our case because the resolution of the synthetic spectrum, used to 
specify the local profiles, is higher than that of the observed spectrum, so that neighbouring 



values o: 



the solution vector P TOt k are expected to be correlated (see Sect. 14.11) . Moreover, 



Barnes! (120041) noted that when the number of basis vectors adopted for the SVD solution is 



sufficiently high, high-frequency oscillations appear in the LSD profile owing to the fitting of 
residual noise. We performed a convolution of our LSD profiles with a Gaussian of FWHM 
of 4.0 km s _1 to overcome this problem, thus obtaining SVD2 profiles in Barnes' terminol- 
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Table 2. Wavelength intervals excluded from the Least Square Deconvolution. 



from 


to 


H 


(nm) 


440.00 


505.00 


514.50 


519.00 


555.00 


557.15 


574.00 


578.20 


588.20 


598.00 


612.50 


632.00 


634.70 


635.70 


654.90 


657.80 


686.50 


790.00 



ogy. Given that our template spectrum has a resolution of 42,000, this treatment does not 
significantly degrade our map resolution. 

Some wavelength intervals have been excluded from the SVD solution because they 
include strong lines or lines affected by the presence of the chromosphere. They are listed in 
Table [2J 

The gain in the signal-to-noise ratio of the extracted mean broadened profile scales 
approximately as the square root of the number of spectral lines included into the LSD 
deconvolution. Considering about 600 — 800 lines in the present application, the gain is 
~ 25, that is we reach a final S/N ~ 2000, on our best exposed spectra. 

3.2 Doppler imaging technique 

The application of Doppler imaging to reconstruct a spot map from a sequence of mean 
broadened profiles, well distributed in phase, requires local mean line profiles for the unper- 
turbed photosphere and the spotted photosphere, respectively We adopted the LSD profile 
of the K5 dwarf Gl 673 (T e g ~ 4500 K) as the unperturbed profile and assumed that the 
profile of the spotted photosphere is a simply scaled version of that of the unperturbed 
photosphere, as discussed in Sect. 13.11 The spectral resolution of the mean local profile is 
42, 000. Note also that Gl 673 has a v sinz negligible in comparison to that of LO Peg, thus 
its mean profile can be considered unaffected by rotational broadening for our purposes. The 
widths of the instrumental profiles of SARG and Elodie are significantly smaller than the 
intrinsic widths of our mean profiles, so the fact that they have been acquired with different 
spectrographs has a negligible effect on our map reconstruction. 

Synthetic mean profiles of LO Peg were constructed for a trial distribution of the spot 
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covering factor / s on the surface of the star, according to the principles outlined in Sect. 13.11 
where / s is the ratio of the area of the spots in a given surface element to that of the given 
element. To construct our maps, we adopted surface elements of size 9° x 9° to oversample 
the actual resolution elements that have an average size of ~ 18°. The parameters assumed 
to synthetize the mean rotational broadened profiles are: a) inclination of the rotation axis 
along the line of sight: % = 45°; b) vsini = 65.84 km s _1 ; c) linear limb-darkening coef- 
ficient u = 0.8; and d) ratio of the intensity in the continuum of the spot mean profile 
to that of th e unsp otted mean profile C s = 0.2. Inclination and v sinz were adopted from 



Barnes et al. 



(120051 ). We show in Sect. 14.21 that small variations of the values of i and v sini 
around the adopted values give systematic residuals in the best fits of th e line profiles thus 

(120051 ). Note that uncer- 



Barnes et al. 



confirming the accuracy of the parameters found by 
tainties in the limb-darkening coefficient and in the detailed shape of the local line profile 
do not significantly affect the reconstruction as they can be almost perfectly compensated 
by minute variations of the vs ini parameter at the level of 0.3 — 0.6 km s _1 (cf., e.g., 



Unruh fc Collier Cameron 



19951 . for details). A variation of the contrast coefficient C s cor- 



responds to a change of the effective temperature of the spots and has an obvious impact 
on the covering factor of the spots. The adopted value corresponds to a spot temperature 
of T efr ~ 3500 K. 

Since LO Peg is a fast rotator, we included the effects of gravity darkening in the con- 
tinuum adopting a stellar radius of 0.78 R , estimated from the vsini, the inclination and 
the rotation period, and a mass of 0.8 M . The effective gravity was computed by adopting 
a simple Roche geometry for the distorted star. 

As it is customary in Doppler imaging reconstruction, the trial spot distribution was var- 
ied to optimize an objective function that consists of a linear combination of the x 2 i obtained 
by fitting the sequence of the LSD profiles, and a regularization function. Here we consider 
Maximum Entropy regu larization (hereinafter ME) according to the approach introduced by 
Collier Cameron! (J1992J). A Lagrangian multiplier specifies the relative weight of the x 2 an d 
of the ME regularization in the objective function. It is determined a posteriori by select- 
ing the maximum value that leads to a distribution of the residuals between the synthetic 
and the LSD profiles that does not deviate by more that one standard deviation from the 
Gaussian distribution ob tained without any regularization. The algorithm of optimization 
is the same used by, e.g.. lLanza et al.l (120021 ) . We extensively tested our code by simulating 
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sequences of line profiles with added Gaussian noise for different spot configurations and 
reconstructing them through the procedure described above. 



4 RESULTS 
4.1 LSD profiles 

An example of the LSD profiles obtained with the procedure described in Sect. 13 .11 is plotted 
in the upper panel of Fig. [U The dotted line is the LSD profile as obtained from the SVD 
solution, i.e., with small ripples due to residual noise fitting, while the solid li ne is the SVD 2 



profile obtained after a convolution with a Gaussian of FWHM of 4.0 km s -1 (jBarnesll2004l ). 
The absolute value of the covariance between the flux value of one of the velocity bin along 
the profile (marked by a vertically dotted line) and those of its neighbour bins is plotted in 
the lower panel. The covariance is obtained with the SVD method before convolving with 
the Gaussian to obtain the SVD2 profile. Note that it is different from zero only in a narrow 
interval of amplitude ~ 6 km s _1 , centred around the given bin, as estimated from the 
FWHM of the covariance curve. Similar values of the FWHM of the covariance curve are 
found also for the other velocity bins along the profile, so we conclude that the actual radial 
velocity resolution of our profiles is ~ 6 km s -1 . The maximum of the covariance curve, 
reached at the radial velocity value corresponding to the given bin, provides us with an 
estimate of the variance of the flux in that bin of the LSD profile. From the value in Fig. 
we find a signal-to-noise ratio of ~ 2000 that confirms our a priori estimate in Sect. 13.11 
The signal-to-noise ratio is further increased by a factor of ~ 2 by the convolution with the 
Gaussian applied to obtain the SVD2 profiles. 

The time series of the residual LSD profiles, with rotation phase plotted versus radial 
velocity, are shown in Fig. [2j We put together all the profiles from the four nights to have 
the best possible phase coverage. Our spectra were taken with an exposure of ~ 1800 s that 
corresponds to their phase width on the plot. The bumps corresponding to three individual 
starspots can be easily traced along the sequence of the profiles. To help follow their motion, 
three curves have been superposed to the plot, each corresponding to one of the spots, 
respectively (see Table |3]). The phase at which a given bump crosses the central meridian of 
the star's disk depends on the longitude of the corresponding spot, while the amplitude of 
its radial velocity oscillation depends on the spot latitude. The fact that the radial velocity 
excursion is limited to approximately ±(45 — 50) km s" 1 and the low inclination of the 
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Figure 1. Upper panel: Example of an LSD profile from a spectrum taken at rotation phase <j> = 0.502 on 13 September 2003. 
The flux, normalized at the continuum level, is plotted versus the radial velocity, measured from the centre of the LSD line 
profile. The width of the radial velocity bins is 2.5 km s . The dotted line is the profile obtained from the SVD solution, 
while the solid line is the corresponding SVD2 profile (see the text). Lower panel: The covariance between the flux of the radial 
velocity bin marked by the vertically dotted line and those of the adiacent bins along the LSD profile at the values of the radial 
velocity reported on the abscissa as computed with the SVD method before applying convolution with the Gaussian. 



Table 3. Latitude ip and phase <j> of the transit across the central meridian of the stellar disk for the three spots corresponding 
to the traces in Fig. [2] The third column lists the linestyle and the colour of the corresponding trace in the same figure. 



if <j> linestyle (color) 

(deg) 



40 0.28 dot-dashed (white) 
50 0.95 dashed (pink) 
65 0.70 dotted (light blue) 



rotation axis of LO Peg (i = 45°) indicate that those starspots are located at high latitudes 
O 40°). 



4.2 Doppler image of LO Peg 

The sequence of the best fits of the LSD profiles obtained with the ME regularization is 
shown in Figs. [31 SI El El and [3 Only 38 profiles were actually fitted since 6 profiles coming 
from low-quality spectra were discarded. The fits are almost always very good, indicating 
that the spectral features associated with the mapped starspots are coherently traced as 
they move across the mean line profiles. 

To show the effect of a small variation of the inclination i and v sin % on the ME best 



fits to the line profiles, we plot in Figs. M and [9] the mean residua 



line profiles for different 



values of those parameters. The values given by lBarnes et al.l (120051 ) are those giving the best 
average fit w hile different values produce larger residuals in some radial velocity intervals 
(cf. Fig. 3 of Barnes et alJboOa ). 
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-75 -35 +35 +75 




2003 September 12-13-14-15 

Figure 2. Time series of the residual LSD line profiles, with rotation phase plotted versus radial velocity. The residuals are 
obtained by subtracting the grand average profile from each profile. White features correspond to starspot signatures (i.e., 
profile bumps). Three radial velocity curves have been traced by eye to help follow the signatures of the most prominent 
starspots. The spot parameters and the linestyles corresponding to each of them are listed in Table [3] respectively. 
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In order to improve the phase coverage, the LSD profiles of all the four nights from 12 to 

15 September were combined together to derive a single ME-regularized map of the surface 
of LO Peg, shown in Fig. [TOj The data are fitted at a reduced y 2 level of 1.1, for a mean 
standard deviation of a = 2.1 x 10~ 4 of the continuum flux obtained for the SVD2 profiles. 
The longitude resolution of the map is limited by the phase blurring due to the exposure 
time and is ~ 18°. The latitude resolution is a function of the latitude itself, with an average 
value of ~ 20° at a latitude of ~ 50°, as derived from our radial velocity resolution of ~ 6 
km s" 1 (cf. Sect. ED). 



Our ME map is similar to those of 



Barnes et al. 



(120051 ). showing a concentration of 



starspots at latitudes higher than 60° with a covering factor varying with longitude. Ap- 
pendages are also found, extending down to a latitude of 10° — 20° and corresponding to 
the moving bumps we traced on the line profile sequence in Fig. [2j Those features are the 
most reliable and their reconstructions are largely insensitive to uncertainties in the mapping 
parameters (cf. Sect. 13.21) because t hey can be directly traced onto the line profile sequence. 



In comparison with the maps of 



Barnes et al. 



(120051 ). our map indicates a slightly higher 



level of mid- and low-latitude activity. However, this result should be taken with caution 
because the latitude resolution of the Doppler imaging technique is lower close to the equator. 
Moreover, ME regularization sometimes tends to produce small concentrations of spots to 
improve the fit of the line profiles, especiall y at phases where observation s are missing (the 
so-called superresolution effect, see, e. g., iNarayan fc Nityanandal Il986l ). The small spot 
just at the threshold of visibility close to the equator at a longitude of ~ 170° may therefore 
be an artefact as suggested by the lack of an obvious corresponding bump in the sequence 
of line profiles and the fact that it quickly disappears when the Lagragian parameter is 
increased. 



5 DISCUSSION AND CONCLUSIONS 



We used the LSD approach of iBarnesI (12004 ). base d on a synthetic spectrum i nstea d of a 



Donati fc Collier Cameron! (119971 ). Our 



line list, as adopted in the original formulation of 
approach takes into account the finite width of the spectral lines due to, e.g., the thermal 
Doppler, pressure and microturbulence broadenings, allowing us a very good fit of the ob- 
served spectra. The advantage is especially apparent in the match of our LSD profiles to 
the continuum in the red and blue extreme line wings. Moreover, the use of the covariance 
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Figure 3. A sequence of LSD profiles (black open diamonds) with superposed their best fits computed with the ME regular- 
ization (red solid line), respectively, fn each panel, the flux normalized to the continuum is plotted versus the radial velocity 
measured from the centre of the unperturbed profile. The rotation phase corresponding to each profile is reported in each panel 
and ranges from 0.0219 to 0.2823 for this subset of profiles. The dimensions of the plotting symbols correspond to one standard 
deviation of the normalized relative flux. 

matrix obtained through the SVD technique, allows us an a posteriori check of the signal- 
to-noise ratio along the average profiles and a determination of their actual radial velocity 
resolution. 

Our Doppler images of LO Peg are similar to those of 



Barnes et al. 



(120051 ). with a polar 



cap and some starspot appendages extending down to « 10° latitude. The spot concentra- 
tions at high latitudes are largely independent of our Doppler imaging assumptions, as the 
moving features associated with them are clearly visible on the time sequence of residual 
line profiles in Fig. [2j The polar cap is not symmetric around the pole, being more extended 
around phase 0.1 and less extend ed around phase 0.5 5. 



Our maps, as well as those of 



Barnes et al. 



(120051 ) . do not show starspots southern than 



~ 10°. However, the presence of small spots or of a band of uniformly distributed spots at 
low latitudes cannot be ruled out on the basis of our Doppler image. The former may escape 
detection due to foreshortening effects, since the inclination of the stellar rotation axis is 
quite low. The latter does not produce any migrating feature on the average line profile, so 



16 N. Piluso et al. 



.000 



0.995 



0.990 - 



.000 



0.995 



0.990 - 




-35 35 -35 35 -35 35 

Radial Velocity (Km/s) 



-35 35 



Figure 4. As Fig. |3] but for phases between 0.3363 to 0.4550. 



any systematic error in the local line profile or stellar parameters can effectively prevent its 
detection. 

The long-term persistence of the polar spot on LO Peg is particularly intriguing when 
its be haviour is compare d with those of other rapidly rotating late- type stars (see Sect. 



4.1 of 



Barnes et al. 



20051 ). The possibility that the polar spot is a non-persistent feature is 



supported by the Doppler maps of AB Dor ( IKiirster et al.lll994l ) that are more numerous and 



better distributed in time than those of LO Peg. The predominance of the Coriolis force on 
magnetic buoyancy and magnetic tension forces on flux tubes emerging from the base of the 
convection zone or the overshoot region at the interface with the radiative core, can explain 
the appearance of flux at intermediate and high latitudes in rapidly rotating stars (see, 



e.g, 



Schiissler et al. 



flows, such mode 



Schrijver fc Title 



19961 ) . In combination with surface diffusion and poleward meridional 
s can explain the pr esence of polar spots and their long lifetimes (e.g., 



2001 



Isik et al.l 



20071 ). The concentration of the magnetic flux toward the 



poles has relevant effects also on theangu 
their magnetized winds (see, e.g., iBuzasi 



ar momentum and mas s loss of such stars through 



1997 



Holzwarth 



20051 ). 
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Figure 5. As Fig. |3] but for phases between 0.4963 to 0.6037. 
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The presence of low-latitude spots is difficult to explain on the basis of a model consid- 
ering slender flux-tubes originating in the overshoot layer or near the base of the convection 
zone and points toward the existence of a distributed turbulent dyna mo working thro ugh- 
out most of the co nvection zone, as in the models considered by, e.g., iBushbyl (120031 ) and 



Covas et al. 



(120051 ). Further support for the presence of a mean-field distributed dynamo 



i n the conv ection z ones of late - type rapidly rotating stars comes from the works of 



Donati 



(1999 


) and 


Lanza 


(2006. 


2007) 
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